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^^ , Context. Only a handful of debris disks have been imaged up to now. Due to the need for high dynamic range and high angular resolution, 
,.1^ • very little is known about the inner planetary region, where small amounts of warm dust are expected to be found, 
tin Aims. We investigate the close neighbourhood of Vega with the help of infrared stellar interferometry and estimate the integrated A'-band flux 
<»>>, originating from the central 8 AU of the debris disk. 

5_j Methods. We performed precise visibility measurements at both short (~30 m) and long (~150 m) baselines with the FLUOR beam-combiner 
"o^ \ installed at the CHARA Array (Mt Wilson, California) in order to separately resolve the emissions from the extended debris disk (short 
C^ baselines) and from the stellar photosphere (long baselines). 

kJ Results. After revising Vega's A^-band angular diameter (0ud = 3.202 ± 0.005 mas), we show that a significant deficit in squared visibility 
. ,-H , (AV^ = 1.88 ± 0.34%) is detected at short baselines with respect to the best-fit uniform disk stellar model. This deficit can be either attributed 
ji^ ■ to the presence of a low-mass stellar companion around Vega, or as the signature of the thermal and scattered emissions from the debris disk. 
^ ' We show that the presence of a close companion is highly unlikely, as well as other possible perturbations (stellar morphology, calibration), 
and deduce that we have most probably detected the presence of dust in the close neighbourhood of Vega. The resulting flux ratio between 
the stellar photosphere and the debris disk amounts to 1.29 ± 0.19% within the FLUOR field-of-view (~7.8 AU). Finally, we complement our 
^-band study with archival photometric and interferometric data in order to evaluate the main physical properties of the inner dust disk. The 
inferred properties suggest that the Vega system could be currently undergoing major dynamical perturbations. 
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1. Introduction a main-sequence (MS) star^ photometri c surveys with IRAS 

dPaiardo-Acosta et all [l999[) and ISO (Laureiis et alJ l2002h 



a 



Vega (HD 172167, AOV, 7.76 pc) is probably one of the most ^^^^ ^^^^^ ^^^^ ^^^^^ ^^^^ ^^ ^3 ^^^^^ ^^^^ significant in 

important stars in astrophysics, as it has b een used as , a pho - ^^^^^^ ^^^^^^ .^ ^j^^ 20 - 25 t^m region. 

tome trie standard for more than a century lIHearnsh aw 1996). „. , ., ,„„„ • ■ , ■, ,, 

., , , ^. ^ , ? — , , Since the mid-1980s, great attention has been paid to Vega 
However, with the advent of infrared space-based telescopes, It , , ,, ,., ^, , , • ,r , •, 
, , , . ,. , , , , ^ andotherVega-Iike stars. They have been imaged from the mil- 
was discovered to have a large infrared excess beyond 12 um ,. • , • , , • ., , ,. • 

, , , . r, rr^ ^1 limetric domain down to the visible, revealing circumstellar 

with respect to its expected photospheric flux dAumann et alJ , .... t- • Tr • i 

r— -Ti i,, . •,.,-, , , , . — '■ n dust arranged m various shapes. For instance, Vega IS known 

119841) . This was identified as the thermal emission from a , , , , , , • •, 

. „ ,. , „ , , , , , n^ «TT ,■ to be surrounded by a smooth annular structure similar to 

circumstellar disk of cool dust located at about 85 AU from , , t^ • „ i • • i n in-->, m r j ^ 

,, ^. , . ^ ,. n ■ „ . .the solar Kuip er Belt, containing about 3 x 10 Mm of dust 

Vega. Since this first discovery of a circumstellar dust around , , „ , , ,ll,r.nol 1^ T 1 1Urvr.^1 i ■ i i i 

grains (Hollan d et aIJII998tlSu et aIJI2005), which also shows 
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Table 1. Individual measurements. Columns are: (1,2) date and time of observation; (3, 4) projected baseline length and position 
angle (measured East of North); (5) squared visibility after calibration and error; (6, 9) HD number of caUbrators used prior 
and after the given data point respectively, means that there was no calibrator; (7, 8, 10, 1 1) quantities used for computing the 
correlation matrix as in Eq. (26) of Perrin (2003): cryi are errors on the estimated visibility of the calibrators. 
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Date 


UT 


baseline (m) 


angle (°) 


(xlOO) 


HD, 


a 


'^v,? 


HDfc 


P 


^^l 


2005/05/21 


06:17 


101.60 


-76.85 


204+ 1.14 





0.000 


0.000 


165683 


0.330 


0.870 




07:31 


127.86 


-90.04 


6.1 ±0.25 


176527 


0.050 


0.870 


176527 


0.060 


1.025 




08:20 


141.07 


-97.43 


2.6 ± 0.08 


176527 


0.026 


1.025 


173780 


0.039 


0.896 




08:59 


148.55 


-102.96 


1.3+0.04 


173780 


0.024 


0.896 


173780 


0.017 


0.895 


2005/05/22 


06:05 


98.63 


-75.25 


23.2 ± 0.22 


159501 


0.240 


0.467 


159501 


0.064 


0.624 




06:24 


105.77 


-79.02 


18.2 ± 0.20 


159501 


0.142 


0.467 


159501 


0.101 


0.624 




06:29 


107.70 


-80.00 


16.8 ±0.18 


159501 


0.120 


0.467 


159501 


0.107 


0.624 




06:39 


111.61 


-81.97 


14.4 ±0.15 


159501 


0.082 


0.467 


159501 


0.115 


0.624 




06:49 


115.39 


-83.83 


12.2 ±0.14 


159501 


0.052 


0.467 


159501 


0.117 


0.624 




06:59 


118.79 


-85.51 


104 ±0.12 


159501 


0.030 


0.467 


159501 


0.115 


0.624 




08:18 


141.45 


-97.68 


2.6 ± 0.07 


173780 


0.014 


0.624 


173780 


0.051 


0.897 




08:23 


142.62 


-98.45 


2.4 ± 0.06 


173780 


0.011 


0.624 


173780 


0.051 


0.897 




08:34 


144.75 


-99.93 


2.0 ± 0.06 


173780 


0.005 


0.624 


173780 


0.049 


0.897 


2005/06/13 


05:22 


33.59 


20.55 


84.2 ± 1.42 


168775 


0.543 


0.152 


168775 


0.362 


0.153 




06:15 


33.85 


13.58 


834 ± 0.92 


168775 


0.269 


0.153 


168775 


0.628 


0.153 




06:46 


33.92 


9.33 


84.5 ± 0.73 


168775 


0.419 


0.153 


163770 


0.523 


0.272 




07:14 


33.96 


5.49 


80.8 ± 0.99 


163770 


0.510 


0.272 


163770 


0419 


0.272 




07:43 


33.97 


1.37 


82.8 ± 1.35 


163770 


0.514 


0.272 


163770 


0.438 


0.272 




08:13 


33.97 


-2.95 


84.5 ± 1.19 


163770 


0.833 


0.272 


168775 


0.129 


0.152 




09:37 


33.82 


-14.54 


83.6 ± 0.67 


163770 


0.123 


0.272 


168775 


0.784 


0.152 




10:04 


33.70 


-18.05 


83.9 ± 0.66 


168775 


0.574 


0.152 


176670 


0.330 


0.167 


2005/06/14 


07:58 


33.98 


-1.26 


85.0 ± 0.90 


176670 


0.521 


0.166 


176670 


0.400 


0.166 


2005/06/15 


06:03 


33.83 


14.15 


844 ± 1.16 


176670 


0.458 


0.167 


176670 


0.457 


0.167 




06:39 


33.92 


9.18 


86.5 ± 1.35 


176670 


0.363 


0.167 


176670 


0.575 


0.166 




07:07 


33.96 


5.37 


84.2 ± 1.41 


176670 


0.544 


0.166 


163770 


0.392 


0.272 



most part of these debris disks, which could potentially har- 
bour warm dust (> 300 K) heated by the star as suggested by 
iFaiardo-Acosta et alJ(ll998l) . Such warm dust would have a sig- 
nature in the near- and mid-infrared that only photometric stud- 
ies have attempted to detect until recently. Indeed, Vega's near- 
infrared (K, L, M) flux was shown to be significantly above the 
modelled photospheric level (Mountain et al. 1985), but this 
discrepancy was most likely due to an inadequate photospheric 
model since Vega's flux is consistent with other A-type stars to 
within standard photometric precision of 2 - 5% (Leaaett et al. 
119861) . In the N band, the best constraint on the thermal emis- 
sion from warm dust has been obtained by nulling interfer- 
ometry, with no resolved emission above 2.1% of the level of 
stellar ph otospheric emission at separations larger than 0.8 AU 
JLiu et alji200 4). At longer wavelengths, the recent measure- 
ments obtained with Spitzer in the far-infrared (Su et al. 2005) 
have not allowed for an investigation of the inner part of Vega's 
disk because of the limited resolution (47 AU at the distance of 
Vega) and because hot dust is not expected to contribute signif- 
icantly to the far-infrared flux. 

In this paper, we use infrared stellar interferometry to in- 
vestigate the inner part of Vega's debris disk. S uch an attempt 
had already been made bv ICiardi et al.l ( 1200 11) . who observed 
Vega with the PTI interferometer on a 110 m long baseline in 
dispersed mode. The poor spatial frequency coverage of their 
observations did not allow clear conclusions, although a sim- 



ple model of a star and a uniform dust disk with a 3 - 6% flux 
ratio was proposed to explain the observations. A more thor- 
ough study of Vega-type stars was performed with the VLTI by 
di Foico et al. (2004), using short and long baselines to sepa- 
rately resolve the two components of the system (stellar pho- 
tosphere at long baselines and circumstellar emission at short 
baselines). Unfortunately, the visibility precision and the avail- 
able baselines at the VLTI only allowed upper limits to be in- 
ferred on the flux of the inner disks. In order to better con- 
strain the near-infrared brightness of Ve ga's disk, we have used 
the sa me method at the CHARA Array (Iten Brummelaar et alJ 
I2OO5I) with an optimised set of baselines. 

2. Observations and data reduction 

Interferometric observations were obtained in the infrared K 
band (1.94 - 2.34 //m) wit h FLUOR, the Fiber Linked Unit 
for Optical Recombination dCoude du Foresto et alJl2003l) . us- 
ing the S1-S2 and E2-W2 basehnes of the CHARA Array, 34 
and 156 metres respectively. Observations took place during 
Spring 2005, on May 21st and May 22nd for E2-W2, and be- 
tween June 13th and June 15th for S1-S2 (see Table0. The 
FLUOR field-of-view, limited by the use of single-mode fibers, 
has a Gaussian shape resulting from the overlap integral of the 
turbulent stell ar image with the fundamental mode of the fiber 
(lGuvonl2002h . Under typical seeing conditions, it has a radius 
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Table 2. Calibrators with spectral type, K magnitude, limb- 
darkened dis k (LP) angular diameter in K band (in m illiarcsec) 
and baseline dBorde et alJ2002HMerand et alJ2005l) . 

S. type A" mag LD diam. (mas) Baseline 



HD 159501 


Kim 


3.14 


1.200 ±0.014 


E2-W2 


HD 163770 


Kllla 


1.03 


3.150 ±0.034 


S1-S2 


HD 165683 


KOIII 


2.9 


1.152 ±0.014 


E2-W2 


HD 168775 


K2 Illab 


1.74 


2.280 ± 0.025 


S1-S2 


HD 173780 


K2III 


2.0 


1.950 ±0.021 


E2-W2 


HD 176527 


K2III 


2.04 


1.765 ±0.024 


E2-W2 


HD 176670 


K2.5 III 


1.6 


2.410 ±0.026 


S1-S2 



of 1" (distance at which the coupling efficiency falls to 3% of 
its on-axis value). 

The FLUOR Data Reduction Software 

JCoude du Foresto et alJ Il997t iKervella et aP l2004 was 
used to extract the squared modulus of the coherence factor 
between the two independent apertures. The interferometric 
transfer function of the instrument was estimated by observing 
calibrators before and after each Vega data point. All calibrator 
stars (TableO were c hosen from two catalogues develope d for 
this specific purpose dBorde et al.ll2002t iMerand et alJ l2005). 
Calibrators chosen in this study are all K giants, whereas Vega 
is an AO dwarf. The spectral type difference is properly taken 
into account in the Data Reduction Software, even though it 
has no significant influence on the final result. The efficiency 
of CHARA/FLUOR was consistent between all calibrators 
and stable night after night to around 85%. Data that share 
a calibrator are afifected by a common systematic error due 
to the uncertainty of the a priori angular diameter of this 
calibrator. In order to interpret o ur data properly, we used a 
specific formalism rPerrin 2003) tailored to propagate these 
correlations into the model fitting process. All diameters are 
derived from the visibility data points using a full model of the 
FLUOR instrument including the spectral bandwidth effects 
(Kervellaetal. 2003) . 

3. Data analysis 

3.1. Stellar diameter 

The measurements obtained with the long E2-W2 baseline 
are particularly appropriate for a precise diameter determina- 
tion, because they provide good spatial frequency coverage of 
the end of the first lobe of the visibility curve (see Fig. [Q. 
Previo us interferometric mea s ureme nts obtai ned in the visi- 
ble by iHanburv Brown et aP (1 19741) and .Mozurkewich et alJ 
1I2OO3I) were used to derive uniform disk (UD) diameters 0ud - 
3.08 + 0.07 {A = 440 nm) and 6Iud = 3.15 + 0.03 {A = 800 nm) 
respectively. In the K band , where the limb-darkening effect is 
not as strong. lCiardi et al.l (l2001 ) estimated the UD diameter to 
be 0UD = 3.24 + 0.01 mas. We have fitted a uniform stellar disk 
model to our E2-W2 data, assuming that Vega's photospheric 
intensity I{(p, A) equals the Planck function with an effective 
temperature of 9550 K for all angles (p. The best-fit diameter 
is 9\jD - 3.218 + 0.005 mas for an effective wavelength of 
2.118 //m, which significantly revises the previously obtained 







3uD = 3.218 + 0.005 mas 
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Fig. 1. Fit of a uniform stellar disk model to the E2-W2 data. 
The quality of the fit is satisfactory (reduced x^ of 1 .29), with 
small residuals that do not display any obvious trend except for 
a small underestimation of the actual data for baselines between 
140 and 150 m. 



estimates'. The quality o f the fit is quite go od (xl - 1.29). 
Unlike in the PTI data of Ciardi et al. (2001), we do not see 
any obvious trend in the residuals of the fit, except for three 
points at projected baselines between 140 and 150 m which are 
slightly above the fit (by ~ 1 .5cr). In fact, Fig.l^not only shows a 
significant discrepancy between the CHARA/FLUOR and the 
PTI data, but also between the 1999 and 2000 PTI data. Ou r 
observations do not support the scenario of Ciardi et al. ( 200 ll). 
who proposed a uniform dust ring with a 3 - 6% integrated flux 
relative to the Vega photosphere in K band to account for the 
trend that they observed in the residuals of the fit obtained with 
a simple limb-darkened disk (LD) stellar model. 

Note that fitting an LD stellar model to our data would 
only marginally improve the fit (see Table O, as the shape 
of the first-lobe visibility curve is not very sensitive to limb 
darkening. Moreover, the actual Umb-darkening parameter may 
be significantly larger than standard tabulated values because 
Vega is suspected to be a fast rotating star viewed nearly 
pol e-on and ffie equatorial darkening may b ias ffie limb pro- 
file taulli ver et al.ll994lPeterson et alJ2 004'). Comple mentary 
observations to our data set, obtained by Aufdenber^ et alJ 
(20061) at ~250 m baselines, confirm this fact and lead to an 
accurate estimation of the TiT-band limb profile, which mostly 
affects visibilities beyond the first null and will not be discussed 
here. 



3.2. Visibility deficit at short baselines 

With this precise diameter estimation, we can now have a look 
at the short-baseline data. In fact, these points do not signifi- 
cantly contribute to the UD fit because of the low spatial fre- 
quencies they sample. Including all the data points in the fitting 



' The A^-band diameter proposed bv lCiardi et alj i200ll) was com- 
puted with the assumption of a flat spectrum for the Vega intensity. 
This explains a large part of the discrepancy with our new value. 
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S1-S2 data Comparison with Ciardi et al. (2001) 
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Fig. 2. The data obtained with the S1-S2 baseline (~34 m) are 
displayed as a function of the projected baseline's position an- 
gle together with the best UD fit computed over the whole data 
set (3.217 mas). The data points are significantly below the best 
UD fit, with a mean visibility deficit AV^ - 2%. The addition 
of a uniform diffuse source of emission in the FLUOR field-of- 
view reconciles the best fit with the data (dotted line). Note that 
there is no obvious dependence of the data points with respect 
to position angle, which would be indicative of an asymmetric 
extended emission. 

procedure gives a best-fit diameter Oud = 3.217 ± 0.013 mas, 
but with a poor;i^^ = 3.36. We show the reason for this poor re- 
duced ;^^ in Fig. 12 where the S1-S2 data points are plotted as 
a function of position angle together with the best UD fit (solid 
line). The observations are consistently below the fit, with a 
Ay2 = 1.88 + 0.34%. 

Systematic errors in the estimation of the calibrator diam- 
eters or limb-darkened profiles are possible sources of bias in 
interferometric observations. In order to explain the measured 
visibility deficit in the S1-S2 data, the diameters of the three 
short-baseline calibrators (Table |3 should have been underes- 
timated by 0.26, 0.35 and 0.33 mas respectively, which repre- 
sent about 10 times the estimated error on their diameters. We 
have made sure that such improbable errors were not present in 
our calibration procedure by cross-calibrating the three calibra- 
tors. No significant departure from the expected LD diameters 
was measured, and the calibrated visibilities of Vega do not de- 
pend on the chosen calibrator. Therefore, it appears extremely 
unlikely that the calibration process may have induced the ob- 
served visibility deficit. 

A limb-darkened stellar model for Vega will not recon- 
cile the best-fit stellar model with the S1-S2 data points (see 
Table |5}, because low spatial frequencies are not sensitive to 
limb darkening. One may think of stellar asymmetry as a pos- 
sible reason for the visibility deficit at short baselines, since 
the position angles of the short and long baselines are al- 
most perpendicular (see Table [Q. However, an oblateness ra- 
tio of 1.07 for Vega would be needed to explain the deficit, 
whi ch would strong ly contr adict previ ous interferometric stud- 
ies (Ivan Belle et aL.2001: .Peterson et aL.2004.) . Other stellar 
features such as spots would not explain this deficit either as 



GuD = 3.202 mas 

+ diffuse emission (1.29%) 

Clardl et ol. 200' (best fit) 
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Fig. 3. Comparison of our E2-W2 data (black dots) with the 
observations of Ciardi et al. (2001) obtained at PTI (triangles: 
data acquired in 1999, squares: data acquired in 2000). The 
data are displayed as a function of spatial frequency, taking 
an equivalent wavelength of 2.145 jum for the FLUOR instru- 
ment (computed fo r a flat stellar spectrum as in the study of 
ICiardi et alJ ( 1200 ih . The Icr errors on the PTI data are shown 
at the bottom of the figure for the sake of clarity. 



they can only appear in the second and higher lobes of the vis- 
ibility function. In fact, a natural explanation to the observed 
visibility deficit would be the presence of an extended source 
of emission in the interferometric field-of-view (e.g. disk or 
companion), which would be resolved with the S1-S2 baseline 
(i.e., incoherent emission). 

In order to assess the amount of incoherent emis- 
sion needed to explain the observed visibility deficit, we 
have added a diffuse emission, uniformly distributed in the 
CHARA/FLUOR field-of-view, to our UD stellar model. 
Fitting this new model to the complete data set gives the fol- 
lowing final result: 0ud = 3.202 + 0.005 mas, /T-band flux ra- 
tio = 1.29 + 0.19%, with a significantly decreased ;t'^ - 1.10 
(instead of 3.36). This result is almost independent of the ex- 
tended source morphology, as the spatial frequency coverage 
of our interferometric data is too scarce to constrain its spatial 
distribution. The extended structure, detected with very good 
confidence (almost Icr), would thus have a relative flux contri- 
bution of 1 .29% with respect to the Vega photosphere in K band 
when integrated over the whole field-of-view (7.8 AU in ra- 
dius). Such an excess does not contradict photometric measure- 
ments in the K b and, which have typical accuracies of 2 - 3% 
(lMegessieiill995 ). The result of the fit is displayed in Fig.|2l 
(dotted line) and Fig.|3l(solid line), as well as in Fig.|4]for a re- 
alistic debris disk model (see Sect. 14. 2t . which gives the same 
best-fit parameters. 

4. Discussion 

In this section, we discuss the possible sources of incoherent 
flux around Vega that could account for the observed visibility 
deficit at short baselines. 
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Table 3. Influence of the limb-darkening parameter a on the 
best-fit diameter and the associated reduced ;^'^ using the whole 
data set, assuming a brightness distribution /(//) - fi" with 
fi - cos 6 the cosine o f the azimuth of a surface element of the 
star JHestroffenl 1 997h . The visibility deficit measured at short 
baselines (S1-S2) with respect to the best-fit model is given 
in the last column, showing a weak dependence on the limb- 
darkening model. 





Best-fit 


X^ 


AV' 


a 


6)ld (mas) 


(all data) 


(S1--S2) 


0.0 


3.217 ±0.013 


3.36 


1.88% 


0.1 


3.264 ± 0.013 


3.14 


1.83% 


0.2 


3.310 ±0.012 


2.96 


1.78% 


0.3 


3.356 ±0.012 


2.82 


1.73% 


0.4 


3.402 ±0.011 


2.71 


1.67% 


0.5 


3.447 ±0.011 


2.64 


1.62% 


0.6 


3.491 ±0.011 


2.61 


1.58% 


0.7 


3.536 ±0.012 


2.60 


1.53% 


0.8 


3.579 ±0.012 


2.62 


1.49% 


0.9 


3.623 ±0.012 


2.66 


1.44% 



from the star For instance, the flux should be increased by 10% 
at 100 mas, by 50% at 200 mas and by 3000% at 500 mas from 
the star in order to reproduce the observed visibility deficit. 
Based on a minimum /T-band flux ratio of 1 .29 + 0.19% and a 
K magnitude of 0.02 for Vega (IMegessien 1 9951) . we deduce an 
upper limit of K = 4.74 ±0.17 for a companion. 

4.1.1. Field Star 

Although Vega is known to be surrounded by a number of faint 
objects {V > 9) wit h low proper motion since the beginning 
of the 20th century (Dommanget& Nvs 2002), these objects 
are far enough from Vega (at least 1') so that they do not inter- 
fere with our measurements. In the i nfrared, neither adap tive 
optics studies (Macintosh et al. 2003; Metc hev et al.l2003l) nor 
the 2MASS survey (Cutri et al. 2003) identified any /T < 5 ob- 
ject within 1' of Vega. In fact, the local density of such objects 
is about 5 x 10"^ per arcmin^ according to the 2MASS survey, 
so that the probability to find aK < 5 source within 1" of Vega 
is smaUer than 4.3 x 10"^. 



4.1. Point source 

Because of our sparse sampling of spatial frequencies, a point 
source located in the FLUOR field-of-view could also be the 
origin of the observed visibility deficit. Regardless of the bound 
or unbound character of the companion, there are essentially 
two regimes to be considered when computing the visibility 
of a binary system, depending on whether the fringe packet 
associated with the companion falls into the FLUOR obser- 
vation window or not. The observation window is defined as 
the total optical path Lqpd scanned by the FLUOR dither mir- 
ror, which is used to temporally record the fringes. The sec- 
ondary fringe packet lies outside the observation window if 
|B a cos 6*1 > Lopd/2, where B is the baseline length, a the 
angular separation of the binary system, the angle between 
the baseline and the orientation of the binary system, and 
^OPD = 102 fira. In that case, e.g. for an angular separation 
larger than 350 mas at a baseline of 34 m, the flux from the 
secondary will contribute incoherently and will lead to the 
same signature as a diffuse emission in the FLUOR field-of- 
view. A binary star with a separation ranging between 350 and 
1000 mas could therefore reproduce the observed visibilities. 
On the other hand, if the secondary fringe packet is inside the 
observation window, it will lead either to a visibility modula- 
tion of twice the flux ratio as a function of baseline azimuth 
if the fringe packets are superposed, or to an enhancement of 
the measured visibility if the fringe packets are separated. Even 
if such behaviour does not seem compatible with the observed 
visibilities, our sparse data cannot definitely rule out a solution 
with a close companion. 

The presence of a point source located within the FLUOR 
field-of-view could thus possibly explain our observations. The 
minimum A'-band flux ratio between the point source and Vega 
is 1.29 + 0.19%, valid for a very close companion (< 50 mas). 
Because of the Gaussian shape of the off-axis transmission, the 
companion would have a larger flux if located farther away 



4.1.2. Physical companion 

At the distance of Vega, the putative companion would have a 
maximum absolute magnitude M^ - 5.15 + 0.17. Assuming 
this companion to be a star of the same age as Vega itself, com- 
prised between 267 and 383 Myr ( Song et al. 2001 ), we use the 
evolutionary models developed bv .Baraffe et aL C1998.) to de- 
duce the range of effective temperature and mass for the com- 
panion: Teff = 3890 + 70 K and M = 0.60 ± 0.025 Mo. This 
roughly corresponds to a n MOV star (IDelfosse et al.-2000.) . 

With a y - /: of 3.65 JBessel & BretJl988h . flie MOV com- 
panion would have a V magnitude of 8.41 and would therefore 
have remained undetected in high resolution visible spectra of 
Vega (M. Gerbaldi, personal communication). Adaptive optics 
studies in the near-infrared would not have noticed the com- 
panion either, due to its very small angular distance from the 
bright Vega (< 1"). At longer wavelengths, the expected in- 
frared excess due to an MOV companion is not large enough 
to be detected by classical photometry as its does not exceed 
2% between 10 and 100 //m. Indirect methods are in fact much 
more appropriate to detect this kind of companion. 

Astrometric measurements of Vega with Hipparcos did not 
detect the presence o f any companio n, with an astrometric 
precision of 0.5 mas (IPerrvmanlll997l) . With a mass ratio of 
4.2 between Vega (2.5AI0) and its putative MOV companion 
(O.6AI0), a 3cr astrometric stability of 1 .5 mas implies that the 
orbital semi-major axis of the putative companion cannot be 
larger than 6.3 mas^ (= 0.05 AU = 4R^,) with a 99% confi- 
dence assuming a circular orbit, which is anticipated for such 
a small separation. Such a close companion, which could also 
fit the interferometric data, would have an appreciable signa- 
ture in radial velocity measurements, unless the binary system 
is seen almost exactly pole-on. Precise measurements recently 



- The astrometric signature of a low-mass companion is given by 
the ratio bet ween the orbital semi-major axis and the mass ratio 
<Perrvmaj2"000) . 
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obtained with the ELODIE spectrometer have shown a relative 
stability ofVega's radial velocity over several months, with am- 
plitudes lower than 100 m/s and a precision of order of 30 m/s 
each (F. Galland, private communication). Assuming that the 
orbital plane of the MOV companion is perpendicular to Vega's 
rotation axis, iiiclined by 5.1° with respect to the line-of-sight 
( Gulliver et alJ 1 1 994) . the companion should be farther than 
80 AU from Vega to be compatible with the measured radial 
velocity stability. In fact, for an MOV companion at 0.05 AU 
from Vega not to display any radial velocity signature at the 
100 m/s level, its orbital inclination needs to coincid e with the 
plane of the sky to within +0.13° ( Perrvmanll2000 ). Even if 
such an inclination is possible, the probability for the system 
to be so close to pole-on is very low (it ranges between about 
6 X 10 "^ and 10"* depending on the assumptions on the sta- 
tistical distribution of low-mass companion orbital planes). In 
conclusion, even though the presence of an MOV companion 
close to Vega could explain the interferometric data, there is 
strong evidence that such a companion does not exist. 

4.2. Circumstellar material 

Circumstellar disks around MS stars are understood to be com- 
posed of second-generation dust grains originating from col- 
lisions between s mall bodies (asteroids) or f rom the evap- 
oration of comets ( iBackman & Parescelll993l) . They are as- 
sumed to be continuously replenished since dust grains have 
a limited lifetime (< 10 Myr) due to radiation pressure, 
Poynting-Robertsori (P-R ) drag and collisional destruction 
JDominik & Decinll2003 o. Several studies have shown Vega 
to harbour a cold ci r cumstellar dust ring ~85 AU in radius 
JHoUand et all Il998t iHeinrichsen et all Il998t iKoerner et all 
l200l[ IWilner et alJl2002h . ISu et al.1 ll2005l) interpreted the ex- 
tended dust emission (up to 600 AU, i.e., 77") detected by 
Spitzer as the signature of dust grains being expelled by radia- 
tion pressure from the Vega system as a result of a recent col- 
lision in the main planetesimal ring and subsequent collisional 
cascade. Even if the presence of dust in the inner part of the 
disk has not been detected yet due to instrumental limitations, 
an equivalent to the solar zodiacal cloud is expected to be found 
around Vega. The thermal and scattered emissions from warm 
grains surrounding Vega could thus be a natural explanation to 
the visibility deficit observed at short baselines, provided that a 
sufficient quantity of dust is present within 8 AU from the star. 
In order to assess the adequacy of a circumstellar disk to re- 
produce the observations, we have fitted our full data set with 
the only known mod el for an inrier debr is disk, i.e., the zo- 
diacal disk model of Kelsall et al.l (ll998V, assuming that the 
inner dust distribution around Vega follows the same density 
and temperature power-laws as for the solar zodiacal cloud. 
The result is displayed in Fig. |4] wherein all data points are 
nicely spread around the best-fit model (as expected, because 
our interferometric data are not sensitive to the particular mor- 
phology of the incoherent emission). The long-baseline data 
are also better fitted than with a simple UD model, because the 
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Fig. 4. Fit of a uniform stellar disk + circumstellar disk model 
to our full data set, using the model of .Kelsall et aL C1998.) . 



presence of the dust disk has som e influence on the slo pe of the 
visibility curve at long baselines ddi Folco et alJ2004l) . The re- 
sulting flux ratio between the whole circumstellar disk and the 
stellar photosphere (1.29 + 0.19%) is the same as with a simple 
model of uniform diffuse emission (Se ct. 13. 2>. with the sam e 
reduced ;t'^ of 1.10. Using the model of Kel sall et alJ (!l998). a 
flux ratio of 1.29% in K band would suggest that the dust den- 
sity level in the inner Vega system is about 3000 times larger 
than in the solar zodiacal cloud. However, we will see later on 
that this model is not appropriate to represent Vega's inner disk 
(it would largely overestimate its mid-infrared flux), so that the 
comparison is not actually pertinent. 

4.2.1 . Physical properties of tine dust grains 

Let us now try to evaluate the main physical properties of the 
dust grains in the inner debris disk. Table |3 gives the photo- 
metric constraints on the near- and mid-infrared excess flux 
around Vega currently available in the literature. Photometric 
constraints at wavelengths longer than 12 jum are not appro- 
priate for our purpose as they are mostly sensitive to the 
cold outer disk (the inner disk is not supposed to produce a 
significant photometric contribution in the far-infrared). The 
large error bars on the photometric measurements take into 
account both the actual error on photometric measurements 
and the estimated accuracy of photospheric models for Vega, 
to which the measurements are compared. Our study is com- 
patible with previous near-infrared measurements but provides 
a much stronger constraint on the inner disk, because inter- 
ferometry spatially resolves the disk from the stellar photo- 
sphere and focuses on the inner part of the disk thanks to the 
small field-of-view. Nulling interferometry at the MMT with 
the BLINC instrument also provides a valuable constraint on 
the mid-infrared excess ( Liu et al.-2004.) . The sinusoidal trans- 
mission map of this nulling interferometer restricts however the 
observation to the part of the disk located farther than about 
125 mas (~ 1 AU) from the star This explains why the result 
of this study is significantly below the estimated mid-infrared 
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Table 4. Available constraints on t he near- and mid-infr ared 
excess around Vega. Re fere nces: (1) Camoins et al. ( 198 a); (2) 
[Blackwell et al.' 0983); (3) 'Rieke et al. (1985); (4) iLiu et alJ 
(2004); (5) Cohen et al. (1992), with th e absolute photometric 
error estimated bv lAumann et al.l (Il984l) . The photometric data 
in references (1), (2) and (3) have been compared to the most 
recent Kurucz photospheric model of Vega (Bohlin & Gilliland 
l2004l) . which has a typical uncertainty of 2% in the infrared 
(this uncertainty has been added to the estimated errors on 
the measurements). Note that the interferometric data from 
FLUOR and BLINC only sample a specific part of the inner 
disk, while the photometric studies include Vega's entire envi- 
ronment. 
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Fig. 5. Map of the x^ as a function of minimum grain size flmin 
and inner radius rn, obtai r ied by fitting the circumstellar disk 
model of lAugereauetalJ (119991) to the SED data of Table |3| 
We have assumed a surface density power-law S(r) oc r^'^ and 
a grain size distribution dn(a) oc a^^-^da, with a maximum size 
of 1500 fim. In this simulation, the disk is composed of 50% 
amorphous carbons and 50% glassy olivines (see text). The 
dashed line represents the distance at which sublimation hap- 
pens for dust grains of a given size (isotherm T - 1700 K). 
The axis labels should therefore read "a in yum" and "sublima- 
tion radius in AU" for this curve. 



photometric excesses, as it is not sensitive to hot grains in the 
innermost part of the disk. 

We have tried to reproduce the Spectral Energy Distribution 
(SED) of the infrared exc ess as listed in T ablelHwith the debris 
disk model developed bv' Augereau et alJ ( 11999 ). For that pur- 

S iose, we took for Vega a NextGen model atmosphere spectrum 
HauschUdt et aDIT999l) with T^s = 9600 K and log(^) = 4.0, 
scaled to match the observed visible magnitude {V - 0.03) 
at a distance of 7.76 pc, which gives a luminosity of 58.7L0. 
Various grain compositions and size distributions were used 
in the disk model, as well as various radial density profiles, 
assuming no azimuthal dependence. In each model, the subli- 
mation temperature of the grains is set to Tjub = 1700 K. At 
a given distance and for a given size distribution, only grains 
large enough to survive the sublimation process can actually 
coexist (see dashed curve in Fig|5}. The normalised differential 
size distribution between a^m and Omax (fixed) is thus truncated 
at Usub, which depends on the radial distance to the star For 
each model, a x^ map is computed for all possible values of 
Qmin (minimum grain size) and ro (inner radius where the disk 
is artificially truncated), adjusting the surface density at ro by 
a least-squares method (see Fig.|5}. The most constraining ob- 
servations in this process are the two interferometric measure- 
ments at 2.12 and 10.6 /im, so that the fitting procedure mainly 
boils down to adjusting the near-infrared flux without produc- 
ing a too strong 10.6 fim emission feature. Comparison of ;t^^ 
values allowed us to infer most probable physical properties for 
the inner debris disk. 

- Size distribution: The inner disk seems to be mainly com- 
posed of hot (-1500 K) and small (< 1 fim) dust grains, 
which emit mostly in the near-infrared. Although larger 
grains (> 10 yum) cannot be ruled out as the main source 



of the excess, such grains generally produce too large a 
mid-infrared flux as they emit more efficiently in this wave- 
length range. This suggests a steep size distribution with 
a small minimum grain size (amin ^ 0.3 yum, assuming 
compact grains). For instance , we find that a size distribu- 
tion similar to that inferred bv lHanneii(ll984l) for cometary 
grains provides a good fit to the SED, a s well as the inter- 
stellar size distribution of Mathis et al. 819771) . Both have 
power-law exponents of -3.5 or steeper On the contrary, 
the size distribution o nGriin et al.l ( ll98S) for interplanetary 
dust particles does not provide a good reproduction of the 
disk's SED, so that the grain size distribution is most proba- 
bly different from that of the solar zodiacal cloud described 
by Reach etal. (2003). 

- Composition: La rge amounts of highly refractive grains, 
such as gra phites (iLaor & Draine.1993.) or amorphous car- 
bons (Zub ko et alJ 119961) . are most probably present in 
the inner disk. This is required in order to explain the 
lack of significant silicate emission features around 10 yum 
(Gaidos & Koresko 2004), which are especially prominent 
for small grains. Silicate grains can still be present in the 
disk, but with a maximum volume rati o of ~70% , using 
the astronomical silicates of Weingartner&Draingj200lk 
or the glassy olivines (Mg2yFe2_2vSi04) o nDorschner et alJ 
( I1995.) with y - 0.5. This is another difference from the so- 
lar zodiacal cloud , which is thought to contain about 90% 
of silicate grains (iReach et al.1 '2003). Such a mixing ratio 
would only be possible around Vega if the grains were suf- 
ficiently big (flmin > 10 yum), SO that the silicate emission 
feature around 10 yum would not be too prominent. 

- Density profile: The inner radius ro of the dusty disk is 
estimated to be between 0.17 and 0.3 AU. Assuming a sub- 
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Fig. 6 . A possible fit of our debris disk model JAugereau et alJ 
1 19991) to the photometric and interferometric constraints of 
Table|4] the diamonds correspond to references (1) and (2), the 
filled circle to this study, the triangle to (3), the square to (4) 
and the cross to (5). The model used here has a size distribu- 
tion dn{a) cc a^^''da with limiting grain sizes a^i^ = 0.1 yum 
and flmax - 1500 yum, a surface density power-law S(r) oc r^'^ 
with an inner radius ro = 0.2 AU, and assumes a disk com- 
posed of 50% amorphous carbon and 50% glassy olivine. The 
solid and dotted lines represent the total emission from the disk 
on a 8 AU field-of-view, respectively without and with the spa- 
tial filtering of interferometric studies, while the dashed line 
takes only the thermal emission into account. The photospheric 
SED, simulated by a NextGen model atmosphere (see text), is 
represented as a dashed-dotted Une for comparison. 

limation temperature of 1700 K, dust grains larger than 
0.5 yum would survive at such distances (see dashed curve 
in Fig. 13 while smaller grains, which are hotter, sublimate 
farther from the star (e.g. at ~0.6 AU for a 0. 1 //m grain). A 
steep power-law for the radial surface density distribution 
has also been inferred from our investigations. A power- 
law exponent of -4 or steeper provides a good fit to the 
SED, as it reduces the amount of dust in the regions far- 
ther than 1 AU and thereb y explains the no n-detection with 
MMT/BLINC reported by'Liu et al.'('2004'). In contrast, the 
zodiacal disk model of Kelsall et al. ( 1998) has a flat sur- 
face density power-law with an exponent around -0.34. 

Using these most probable parameters for the inner 
disk and a mixed composition of 50% amorphous carbons 
llZubkoetal.lll99d ) and 50% glassy olivines (MgFeSi04, 
IPorschner et alJ ll995). we have obtained a relatively good fit 
to the SED as illustrated in Fig.|6l where we see that the ther- 
mal emission from the hot grains supersedes the contribution 
from scattered light at wavelengths longer than 1.3 /vm. Based 
on our model and assuming a size distribution dn{a) oc a^^'^da 
with (Jmin = 0.1 yum and cimax - 1500 fim, we can deduce 
estimations for the dust mass in the inner 10 AU of the disk 
(Aldust ~ 8 X lO^^A^e, equivalent to the mass of an asteroid 
about 70 km in diameter) and for the bolometric luminosity ra- 
tio between the inner disk and the star (Ldisk/^* ~ 5 x 10 "*). 



Because of the high temperature of the grains, the luminosity of 
the inner disk is more than one order of ma gnitude larger than 
the luminosity of the outer disk estimated bv lHeinrichsen et alJ 
( 1998), even though it is almost 10^ times less massive than 
the outer disk. These results need to be confirmed by future 
studies, as the SED of the inner disk is still relatively poorly 
constrained. They have been included in this paper to demon- 
strate that the presence of warm circumstellar dust can repro- 
duce the various observations, and to provide a plausible dust- 
production scenario as discussed below. 



4.2.2. A possible scenario for tine presence of iiot dust 

In fact, three main scenarios may explain the presence of small 
dust grains so close to Vega. As in the case of the solar zodia- 
cal cloud, they could be produced locally, e.g. by collisions be- 
tween larger bodies arranged in a structure similar to the solar 
asteroidal belt. Another local source of small grains is the evap- 
oration of comets originating from the reservoir of small bodies 
at ~85 AU from Vega or fro m an inner population of icy bod- 
ies as in the case of/? Pic JBeust & Mo rbidelli 2000). Finally, 
dust grains produced by collisions in the outer disk could drift 
towards the inner region because of P-R drag. However, this 
latter scenario cannot be c onnected to the recent collision(s) in 
the outer disk suggested by'Su et al'flOOy), because of the long 
timescale of P-R drag (2 x 10^ yr. Dent et al. 2000). Moreover, 
due to the much shorter collisional timescale (5x10^ yr in the 
outer disk), this process is not very efficient and is therefore un- 
likely to produce the observed amount of dust in Vega's inner 
system. Our observations cannot discriminate between the two 
remaining scenarios, even though a cometary origin is favoured 
by the steep size distribution of dust grains cHanner.1984.) and 
by the small inner disk radius. 

Due to radiation pressure, small grains will not survive in 
the Vega inner disk m ore than a few year s before being ejected 
toward cooler regions JKrivov et alJ2000 ^. Larger grains would 
survive somewhat longer, but not more than a few tens of years 
due to the high collision rate in the inner disk. A large dust pro- 
duction rate (~ lO^'^AIe/yr) is thus needed to explain our ob- 
servations, suggesting that major dynamical perturbations are 
currently ongoing in the Vega system. An attractive scenario 
would be an equivalent to the Late Heavy Bombardment that 
happened in the solar s ystem in the 700 Myr following the for- 
mation of the planets ( iHartmann et alJ l2000). i.e., at a period 
compatible with the age of Vega (~350 Myr). Such a bombard- 
ment, most probably triggered by the outward migration of gi- 
ant planets JGomes et alJl2005l) . could explain the presence of 
small grains around Vega both in its outer disk, due to an en- 
hanced collision rate in this part of the disk, and in its inner 
disk, due to the high number of comets sent toward the star by 
gravitational interaction with the migrating planets. Although 
the presen ce of giant pl anets around Vega has not been con- 
firmed yet, Wvatt (2003) has suggested that the outward migra- 
tion of a Neptune-sized body from 40 to 65 AU could explain 
the observed clumpy structure in Vega's outer disk. 
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5. Conclusion 

In this paper, we have presented high precision visibiUty mea- 
surements obtained on Vega at the CHARA Array with the 
FLUOR beam-combiner. The presence of a significant deficit 
of visibihty at short baselines with respect to a simple uni- 
form disk stellar model led us to the conclusion that an ad- 
ditional source of /T-band emission is present in the FLUOR 
field-of-view centred around Vega (1" in radius), with an esti- 
mated excess of L29 + 0. 1 9% relative to the photospheric emis- 
sion. Among the possible sources for this excess emission, the 
presence of dust grains in the close vicinity of Vega, heated by 
the star and radiating mostly in the near-infrared, is proposed 
as the most probable one. Vega, a prototypical debris-disk star 
surrounded by a large quantity of dust at about 85 AU, was al- 
ready suspected by several authors to harbour warm dust grains 
arranged in an inner circumstellar disk. Previous studies were 
however limited to a precision of a few percent on the total in- 
frared flux of the Vega system and therefore did not provide a 
precise estimation of the excess emission associated with the 
inner disk. 

Thanks to our precise determination of the integrated K- 
band emission emanating from the inner 8 AU of the Vega 
debris disk, we were able to infer some physical properties 
of the dust, which is suspected to be mainly composed of 
sub-micronic highly refractive grains mainly concentrated in 
the first AU around Vega and heated up to 1700 K. An esti- 
mated dust mass of 8 X 10"^ Ale and a fractional luminosity of 
~ 5 X 10"^ are derived from our best-fit model. We propose 
that a major dynamical event, similar to the solar system Late 
Heavy Bombardment, might be the cause for the presence of 
small dust grains in the inner disk of Vega. 
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